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Abstract. Motivated by recent studies in geophysical and planetary sciences, we investigate the 
PDE-analytical aspects of time-averages for barotropic, inviscid flows on a fast rotating sphere S 2 . 
Of particular interests are the incompressible Euler equations. We prove that the finite-time-average 
of the solution stays close to a subspace of longitude-independent zonal flows. The intial data can 
be arbitrarily far away from this subspace. Meridional variation of the Coriolis parameter underlies 
W) this phenomenon. Our proofs use Riemannian geometric tools, in particular the Hodge Theory. 

j— ( 1. Introduction 

Recent studies have seen increasing understandings of global characteristics of geophysical flows 
on Earth and giant planets in the Solar System. Simulations and observations have persistently 
shown that coherent anisotropy favoring zonal flows appears ubiquitously in planet scale circula- 
tions. For a partial list of computational results, we mention [5] for 3D models, [201 ES El El HI] 
for 2D models, and references therein. These highly resolved, eddy-permitting simulations are 
made possible by rapid developements of high performance computing. On the other hand, we 
have observed zonal flow patterns (bands and jets) on giant planets for hundreds of years, which 
^><0 has attracted considerable interests recently thanks to spacecraft missions and the launch of the 

Hubble Space Telescope (e.g. [7], [E]). Figure 1 shows a composite view of the banded strucure of 
Jovian atmosphere captured by the Cassini spacecraft ([12]). There are also observational data in 
the oceans on Earth showing persistent zonal flow patterns (e.g. [T6l [T71 [TT] ). 

Remark 1.1. In the existing literature, a necessary process for the zonal flow pattern to emerge is 
averaging of data and/or simulations over a period of time usually decades long. Also, a number of 
heuristic arguments (e.g. [6],) are made pointing to the north-south gradient of Coriolis parameter 
as the underlying machanism, even in the absence of temperature/density gradient and vertical 
variability. 

To this end, we study inviscid, barotropic geophysical flows on a unit sphere § 2 centered at the 
origin of R 3 and fast rotating about the z-axis with constant angular velocity. Let vector field 
u(t, q), tangent to S 2 at every point q £ E 2 , denote the fluid velocity relative to this rotating 
frame. Throughout this paper, we represent any point q S S 2 either by its relative-to-the-frame 
cartesian coordinates (x, y, z) or its relative-to-the-frame spherical coordinates (8, <j>) with 9 being 
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Figure 1.1. This true-color simulated view of Jupiter is composed of 4 images taken 
by NASA's Cassini spacecraft on December 7, 2000. Credit: NASA/JPL/University 
of Arizona [12] . 

the colatitude and 4> the longitude^] Also, let be the unit vector in the zonal direction of 
increasing longitude and eg be the unit vector in the meridional direction of increasing colatitude. 

In this study, we focus on a canonical PDE system: the incompressible Euler equations under 
the Coriolis force ([HUE]), 

5(U + V u u + VP = V, divu = (1.1) 

where constant e, called the Rossby number, scales like the frequency of the frame's rotation 
(usually 0.01 ~ 0.1 at a global scale), and u -1 denotes a counterclockwise 7r/2-rotation of u on S 2 . 
Cartesian coordinate z indicates how the Coriolis parameter varies along the meridional direction. 
Note that the Coriolis force is not uniformly large and actually vanishes on the equator. It is the 
large gradient of the Coriolis parameter that drives the zonal flow patterns. 

Note that results have been established concerning solution regularity of the above systems and 
related ones in the fast rotating regime with £<1. Please refer to [21 [31 [8] and references therein 
for further discussion. 



At z = 1, we fix 9 = but (f> can be arbitrary. Such singularity issue does not occur for cartesian coordinates 
a;,y, 2 GC7° c (§ 2 ). 
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Our theoretical investigation is then focused on the fast rotating regime with e -C 1 and in 
particular, the nature of the time- averages of u: 

u(T,.):=^j\(t,.)dt (1.2) 

for positive times T. 

The main result is stated as following. 



Theorem 1.1. Consider the incompressible Euler equations (1.1) on S 2 with div-free initial data 
u(0, •) € H k (E> 2 ) for k > 3. Define the time-averaged flow u as in (1.2). Then, there exist e- 
independent constants Cq,Tq s.t. for any given T G [0, To/||uo||#fc], there exist a function /(•) : 
[—1,1] > K and a universal constant C s.t. 

u(T,x,y,z)-V ± f(z)\\ Hh _ 3 ^ < C e(^ + M 2 ). (1.3) 

Here, Mq := \\uo\\ H k indicates the size of initial data. In spherical coordinates, the approximation 
V x /(z) is 

V ± f(z) = -/'(cos0)sin0e0. 

Our theoretical result proves computational and observational results in the literature mentioned 
at the beginning of this paper, especially Remark |1.1[ Note that our result shows that the zonal- 
flow pattern becomes prominant with decreasing Rossby number e \ 0. In other words, the 
time-averaged flow u is only 0(e) away from a very restricted subspace consisting of longitude- 
independent zonal flows. The initial data, on the other hand, do not need any filtering and can 
be arbitrarily far away from that subspace of zonal flows. Our proofs below will suggest that such 
unique pattern is essentially due to the Coriolis parameter z/e that varies meridionally from the 
strongest at the poles to zero on the equator. 

Rossby number e is typically at magnitude 0.01~0.1 for Earth oceans, which results in a time 
scale of magnitude 10~100 Earth days according to Theorem |1.1[ This suggests that zonal flow 
patterns can occur at time scales far below those used in the literature. In fact, the Rossby number 
is even smaller for giant planets, leading to the direct observability of banded structures. 

Our result for rotating incompressible Euler equations uses the abstract framework of the fol- 
lowing lemma. 

Lemma 1.1. Consider time- dependent equation 

u t = ^£[u] + f(t,q) 

over certain spatial domain 0. Here, e > is a scaling constant, f(t,q) a source term and C a 
linear operator independent of time. Let operator IT denote (some) projection onto the null space 

null{£} 

of C. Assume a priori u, /(t, q), C[u], II [u] have enough regularity as needed. 

null{£} 

Then, under the assumption 

Il u - n n {£} u ll^i(f7) < c\\c[u}\\ Hk2{n) (i.4) 



for some constant C , holds true the following estimate on the time-average of u, 

H k i ~ \ T 
H k 2 and M' : = max t6[QiT ] \\f(t, -)\\ H k 2 . 













[ U udt 


T 




/O null{£} 



where constants M : = max te [ 0) T] || u (^ 



Remark 1.2. The key hypothesis (1.4) is automatically true in a finite- dimensional space if u is a 
vector in W 1 , £ a linear transform W 1 i— > IR n and H the I 2 -projection onto null{£}. In such case, 

null{£} 

hypothesis (1.4), with the norms understood as I 2 norm on both sides, amounts to the boundedness 
ofC- 1 : image{£} R n /null{£} 

Proof of Lemma \l.l\ First, transform the original equation into 

1 



-C[u 



n u] + / 

null{£} 



and apply time-averaging ^ J Q • dt on both sides 

r-T 



1 



(u(T,.)-u(0,-)) 



1 

ef 



n u]dt + - 



fit, -)dt 



Since all necessary regularities were assumed available and C was assumed to be linear and inde- 
pendent of time, we argue that ■ dt and £[■] commute, so that the above equation becomes 



1 



(u(T,-)-u(0,-)) 



1 



udt 



1 

f 



V : 
Due to the factor - in the first term on the RHS, we have 



IT udt 

ull{£} 



+ 



T 



f(t,-)dt. 



C 



udt 



1 
T 



n udt 

null{£} 



< £ 



H k 2 



(2M 



+ M' 



Finally, apply estimate (1.4) to arrive at the conclusion. 



□ 



Note that, the constant M used in the above lemma depends on size of the solution up until 
time T and is not necessarily independent of e. A priori estimates uniform in e are therefore in 
order. The proof requires considerations beyond the well established energy methods, which will 
be explained in Appendix B. 

Proof of Main Theorem \l.l\ Having Lemma 1.1 and e-independent estimates in Theorem |6.1[ it 
suffices to study properties of properly defined operators £ (c.f. Definition 2.1) and IT (c.f. 

null{£} 

Lemma 3.2), and to finally prove estimate (1.4) (c.f. Theorem 4.1). □ 

The rest of the paper is organized as following. We start Section 2 with describing a version of 
the Hodge decomposition in terms of differential operators on E 2 . The definitions of these operators 
are given the Appendix. An elliptic operator £ that plays the same role as the £ in Lemma |1.1| 
is defined by the end of Section 2. In Section 3, we characterize the null space of £, identifying 
null{£} as the space of longitude-independent zonal flows (c.f. Lemma 3.1). We also define the 
projection operator n and its complement. In Section 4, we obtain Sobolev-type estimates, in 

null{£} 



particual (1.4), regarding C and II using the spherical coordinates and spherical harmonics. In 

' ' null{£} 

Appendix A, we give the rigorous definitions of differential operators on surfaces such as S 2 and 
prove related properties. It is necessary to adopt coordinate-independent differential geometric 
tools since any global coordinate system on § 2 is bound to have singularity issues. On the other 
hand, one can formally use spherical coordinates as well as cartesian coordinates for most of the 
arguments presented in this paper, knowing their validity is justified. In Appendix B, we prove an 
e-independent estimates after carefully examining commutability properties of some differential- 
integral operators on a sphere. 



2. Hodge Decomposition 

The Hodge decomposition theorem ([21], |19j ) confirms that for any fe-form u on an oriented 
compact Riemannian manifold, there exist a (k — l)-form a, (k + l)-form /3 and a harmonic A;-form 
7, s.t. 

lj = da + 5(3 + 7. 

In particular, if the manifold is a surface in the cohomology class of S 2 (loosely speaking, there 
is no "hole" or "handle"), then there exist two scalar- valued functions <3? (called potential) and \P 
(called stream function) such that 

u = Ui rr + u; nc where Ui rr := V<3? and Uj nc := V~ L, I / . (2-1) 



Please refer to (5.14) in the Appendix and the discussion that leads to it. 
Moreover, the decomposition satisfies 

CUrl Uj rl . = div Ui nc = 

and 

u^r = VAga 1 div u, u inc = V ± A~2 curl u. (2.2) 

In other words, any (square-integrable) vector field on § 2 can be written as superposition of an 
incompressible and an irrotational vector fields that are determined by ( |2.2[ ). Such decomposition 
is unique because a harmonic scalar-valued function on a sphere (and any surface in the same 
cohomology class) is always constant and therefore A^ 1 is unique up to a constant. 

For simplicity, we will use A for A§2 from here on. Also, we assume that, unless specified 
otherwise, 

A -1 / always has zero global mean over S 2 . (2.3) 

We postpone the differential-geometric definitions and properties of VjV^, curl, div, A on S 2 
till the Appendix. 



Now, rearrange (1.1) as 

z 



u x - V u u = d t u + VP. 

e 
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Observe that on the RHS, c^u is incompressible and VP is irrotational. Thus, the RHS is the 



unique Hodge decomposition of the LHS, which satisfies the elliptic PDEs (2.2). In particular, the 
incompressible part dtu is uniquely determined by 



d t u = V- 1 A -1 curl 



U" 



V u u 



(2.4) 



This is indeed an equivalent formulation of the original incompressible Euler's equation (1.1). 
In the context of Lemma |1.1[ we define the following operator 



Definition 2.1. For any u, not necessarily div-free, define 

C[u] := V ± A~ 1 curl (zu J 



(2.5) 



Here, A 1 follows the convention (2.3). 



Then, (2.4) can be reformulated as, 



d t u + V^A" 1 curl (V u u) = -C[u\. 
3. Null Space of C and associated L 2 -orthogonal projection 



(2.6) 



We observe that the definition (2.5) naturally implies a sufficient condition for C[u] = is 

curl (air 1 ) = 0. 

Indeed this is also necessary by the virtue of ( |5.13 ), namely curlV^ = A, so that 

curl£[u] = curl V^A^ 1 curl (zu 1 ) = curl^u -1 ). 
Thus, for velocity field u, not necessarily div-free, 

£[u] = curl^u- 1 ) = 0. 

Further analysis reveals the following lemma. 
Lemma 3.1. (Characterization of null{£}) For div-free u with sufficient regularity, 

C[u] = u = V ± g(z) = -g'(cos9)sin6ej ) 



(3.1) 



(3.2) 



for some function g : [—1,1] i— >■ R. Thus, we identify null{£}, when restricted to div-free velocity 
fields, with the space of longitude-independent zonal flows. 



Proof of Lemma 3.1. Apply the product rule (5.16) in the Appendix to curl (zu ), 

curl^u- 1 ) = (Vz) • u + zdivu = (Vz) ■ u (3.3) 

where the • product is given by the natural metric on § 2 induced from M 3 and we used the in- 
compressibility condition divu = 0. Therefore, by (3.1) and (3.3), the null space of C is identified 

as 

For any div-free u , £[u] = (Vz) • u = 0. (3.4) 
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Since is in the meridional direction, (3.4) implies any u in the null space of C flows in the 
zonal direction. But there is more than that. Hodge decomposition (2.1), (2.2) implies 



div u = 



u = V ^(x,y,z) 



(3.5) 



with being a unique scalar function with zero global mean. Combining (3.5) it with (3.4), we 
have, for any incompressible velocity field u, 

u = V ± ^ G null{£} (Vz) || (V*). 

The condition (Vz) || (V^) implies ^ is a function of z only. Thus, we arrive at the conclusion. 



The very last term in (3.2) is due to the fact that, in spherical coodinates, 

V^z = — sin^e^. 

□ 

It is then easy to show the following characterization of II , the L 2 -orthogonal-projection 

null{£} 

operator onto null{£}. 

Lemma 3.2. (Characterization of H ) For any div-free vector field u E L 2 (§ 2 ), its L 2 -orthogonal- 

null{£} 

projection onto null{£} satisfies 



n u 

QU11{£} 



u • e r , 



1 



2-7T sin 9 



u • e d 



fC{9) C <P ' c Zjn 0111 u yC{9) J 

where §qiq\ is the line integral along the circle C{9) at a fixed colatitude 9. 



(3.6) 



Several remarks are in order. First, among all possible projection operators, we chose one 
that nullifies the L 2 (S 2 ) orthogonal complement of null{£}. Secondly, intuitively, II u at a 



null{£} 



given latitude is a uniform zonal flow equal to the mean circulation of u at that latitude; thus, 
(id — II )u is of zero circulatioi 

null{£} 

with the orthogonality condition 



(id — II )u is of zero circulation along the circle at a fixed latitude. Such intuition is consistent 

V null{£} 



n u) ■ (id - n W = o. 

2 'nulKC} null{£} 



Thirdly, even though (3.6) runs into singularity at the poles 9 = and 9 = n, such singularity is 



removable. In fact, apply the Stokes' theorem on the RHS of (3.6) so that, for 9 G (0,7r), 

1 



n u 

null{£} 



interior of C(9) 



2tt sin 9 

and by taking the limit as 9 — > 0+ and 9 — > n— , we obtain 

= 0. 



curlu e„ 



n u 

null{£} 



poles 
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In terms of the stream function, for any div-free velocity field u = V _L, I / which amounts to 
u = (dg^)eg 1 - + ^^^^ e^- 1 with removable singularity at the poles, 

jy v± *> - h Of 8s * *) e * = v± (s f * # ) ■ (3 - 7) 

In other words, II maps the stream function to its zonal means. 

null{£} 



4. Key Estimates 



This section is dedictated to proving an estimate similar to (1.4). A convinient tool in studying 
Sobolev norms of functions on § 2 is the spherical harmonics. 

To this end, introduce the spherical harmonc of degree I and order m, 

Yr((f>,6) = N^e^QY 1 (cos 6), for I = 0,1,2,... and m = -I, -I + 1, 0, I - 1,1 



where the normalizing constant 2VJ™ = "J ^^p- so that || Y5 m ||£2(g2) = 1. It satisfies the 

eigenvalue problem 

AY; m = -1(1 + 1)17" (4.1) 

and orthonomal condition 

(ym^Y™ ) L 2( S 2) = 5w5 mm i. 

Here and below, 

(f,g)L 2 (& 2 ) : = / Jgdtt 
Js 2 

with dQ being the area element of §> 2 and locally equals sinOdOdcf). The assoicated Legendre 
polynomial QV l (z) satisfies the general Legendre equation 



dz 

and can be expressed via the Rodrigues's formula, 

flPW->-^)»^(i-^- 

In order to estimate the Sobolve norms (esp. H norms) of Y™, we take the L 2 (§ 2 ) inner product 
of (4.1) with Y (omitting indices for simplicity), invoke Green's identity (5.15) to calculate 

1(1 + 1) =l(l + l)(Y,Y) L2m 

= -{Y,AY) Li§2) = (VY,VY) L 2 (s2) 

which implies 

Y G H 1 ^ 2 ) so inductively Y G H k (§ 2 ), k > 0. 

As a matter of fact, a little more rigor is needed in defining Sobolev norms on a manifold, but we 
will skip the technical details and only use the fact that 

C' k (l + l) k < \\YT\\hH&) < C k (l + l) k , for Z>0, k>0 



This estimate allows us to use the following definition, among many other equivalent definitions 
(|19|), of the H k norm of a scalar- valued function * defined on § 2 . 



Definition 4.1. For a scalar function ^ with J g2 * = and series expansion 

oo I 

* = E E cir, ^re c = (*,Yr) LH ^), 



(4.2) 



i=l m=—l 



define it H k norms, among other equivalent versions, as 



1*1 



H k :- 



oo I 



v E E( 1 + / ) 2fc i^ 
\ i=i 



m\2 



(4.3) 



m=—l 



Remark 4.1. i/ere and below, we always start series and sums with I = 1 and assume i/jq = 
/** = <). 

Consequently, we define H k norms for u. 
Definition 4.2. For a vector field u with Hodge Decomposition 

u = V* + V" 1 * u/itfj J s a * = J g2 * = 0, 
we define its H k norm, among other equivalent versions, as 

; (4-4) 



I u IIh* := 



l$ll 2 , , , + II*" 2 



In particular, if u is div-free with u = V" Ll I / and L 2 * = 0, th 

V x * 



en 



oo I 



m|2 



(4.5) 



m=—l 



Remark 4.2. #ere and below, we always choose *, * tuzi/i zero global mean such that the above 
definition is consistent with \\0\\ H k =0. 

Remark 4.3. Apparently, under above definition, div-free and curl-free vector fields are H k - 
orthogonal. 

We now characterize operator C using the spherical harmonics. Let incompressible velocity field 

u = V- 1 *. 



Then, by definition ( |2.5[ ) and identity ( |3.3[ ) 

£[V X *] = V^A" 1 (Vz • V x *) 
It is easy to verify that, in spherical coordinates, 

Vz = -sin9e e and V" 1 * = (<9 e *)e/ + 
Thus, combining the three equalities above, we obtain 



sm( 



(4.6) 



Lemma 4.1. (Spherical-harmonic representation of C) For a scalar function ^ with a series 
expansion (4.2), the identity (4.6) leads to 



oo I 



^ ^ Z(Z + 1 

2=1 m =-Z v 



(4.7) 



Here, we used the fact that d^Y™ = iraY™ and A 1 Y^ n 



1 vm 



1(1+1) A I 



Y, m for I > 1. We exclude 



/ = from the series due to V'o = J§2 ^ = 0. We also exclude m = since it doesn't contribute to 
(4.7) anyway. 

We now use spherical harmonics to characterize the projection operator II given in (3.7). It 

null{£} ■ ■ 

follows from (4.7) that 



uEnull{£} iff u = V ± X;E^/ m 1 



1=1 m=0 



which is consistent with Lemma 



3.1 



since Y,° is a function of 9 only. 



Thus, the only modes that survive II are those with m = since II is an L 2 -orthogonal 

null{£} null{£} 

projection and 



A ymym' =l(l + l)5 iv 8 mml . 



Lemma 4.2. (Spherical-harmonic representation of II .) For a scalar function ^ with a series 

null{£} 

expansion (4.2), 



OG 



(id- n )(v- L *) =v- 

null{£} 



2=1 m=0 

oo I 



E E 



Z=l m=—l 



\ 



/ 



(4.8) 



(4.9) 



Note that the above 2 equations can also be derived from (3.7) together with the fact that 

r27r 



1 

2^ 



o 



Y l m ((f>,9)d<l> = 8 m o. 



Combining (4.7) with (4.9) and using the absence of m = modes from both series, we induce 
that, when C is restricted to the image of (id — II ), its null space is trivial and its inverse is 

null{£} 

"bounded" (as noted in Remark 1.2, this is automatically true for linear transform C : W 1 i— > W 1 ). 
More precisely, 



Theorem 4.1. For any div-free vector field u G H k (S 2 ) and k>0, 



u- n u 

null{£} 



< 



c\u- n ul 

L null{£} J 
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|£[u] 



Proof. Consider the stream function ^ so that u = V 1 4 r . Combining (4.5) and (4.9), we obtain 



u- n u 

null{£} 



OO I 



E E( 1+ o 2(fc+1) h/> 



m 1 2 



m=—l 



Combining (4.5) and (4.7), we obtain 



\C\u] 



\ H k+2 



OO I 



EE (i+o 

\ z=i m =-J 



2(fc+3) 



1(1 + 1) 



The key observation here is that m = modes are absent in both series; thus, by a simple inequality 

2 



(l + f)2(*+l) < (l + i)2(*+3) 

we arrive at the conclusion! 



m 
1(1 + 1) 



for I > 1, Iml > 1, 



□ 



5. Appendix A: Preparation in Differential Geometry 

Let M denote a 2-dimensional, compact, Riemmanian manifold without boundary, typically the 
unit sphere M = S 13-1 endowed with metric g induced from the embedding Euclidean space M 3 . 
Let p G M denote a point with local coordinates (pi,P2)- Any vector field u in the tagent bundle 
TM is identified with a field of directional differential operator which is written in local coordinates 
as 

5 ° dpi ' 



We use the notation 



V U /:=E< 



10/ 
5pi 



to denote the directional derivative of a scalar- valued function / in the direction of u. Using the 
orthogonal projection Pmi T ^3^xM induced by the Euclidean metric of M 3 , we define the covariant 
derivative of a vector field v £ TM along another vector field u £ TM, 



V u v := Proj TR3 ^ TM E( V u^) e i- 



(5.1) 



i=l 



Here, v is expressed in an orthonomal basis of I 3 as v = ui e l + «2e2 + ^363. 

The metric g is identified with a (0,2) tensor, simply put, an 2 x 2 matrix (5^)2x2 i n local 
coordinates. Thus, the vector inner product follows 

d d 

9(-^-^^—)=9ij forl<z,j<2. 
dpi opj 
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5.1. Hodge Theory. f[2Tl[19]l 

The Hodge *-operator, defined in an orthonormal basi^ satisfies 

*dp\ = dp2, *dp2 = —dpi, *1 = dpi A dp2, *(dpi A dp^) = 1. 

Using the Hodge star operator, we define the co-differential for any £;-forms a in an n-dimensional 
manifold, 

codifferential : 5a := (-l) k *~ 1 d*a = (-l)"^ 1 ^ * d * a, 
and in particular, for n = 2, 

5a = — * d * a. 

In the case when the Riemannian manifold M has no boundary, the codifferential is the adjoint of 
exterior differential w.r.t. L 2 (M) inner product induced by the given metric g, 

(da,f3) L 2 {Mtg) = (a,6P} L 2( M) gy (5.2) 

The Hodge Laplacian (a.k.a. Laplace-Beltrami operator and Laplace-de Rham operator) is then 
defined by 

A H :=d5 + 5d. (5.3) 
In particular, for a scalar- valued function / in a local basis |^"| with metric g, it is identified as 

& H f = --^=Y,d i (^\g~\g i id j f) 
\/\q\ " 

where (g lJ ) is the matrix inverse of (gij). Thus, on a surface M, the Hodge Laplacian Ah defined 
in(5.3) amounts to -1 times the surface Laplacian Am- In particular, if M is a two-dimensional 
surface, then 

for scalar function /, Am/ = —5df = *d* df (5-4) 

since 5f = on a two-dimensional manifold. 
For now on, we will use A for Am- 

The Hodge decomposition theorem in its most general form states that for any /c-form oj 
on an oriented compact Riemannian manifold, there exist & (k — l)-form a, (k + l)-form /3 and a 
harmonic fc-form 7 satisfying A// 7 = 0, s.t. 

oj = da + 5/3 + 7. 

In particular, for any 1-form w on a 2-dimensional manifold with the 1st Betti number (loosely 
speaking, there is no "holes"), there exist two scalar-valued functions \E' such that 

u = d$ + = d$ - (5.5) 

Here, we used the Hodge theory to equate the dimension of the space of harmonic fc-forms on M 
with the fc-th Betti number of M. For the cohomology class containing the unit sphere S" 3-1 , the 
0th, 1st and 2nd Betti numbers are respectively 1, 0, 1. 



2 The existence of such basis is guaranteed by the Gram-Schmidt orthogonalization process. 
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5.2. In Connection With Vector Fields. Let b and (J denote the musical isomorphism between 
TM and T*M induced by the given metric g, i.e. for vector field a = s ^ i a l di and covector field 
(1-form) A = Y^ i Aidp i 

index-lowering, C^^a % dif = gjja l dpj 

index-raising, Aidp 1 )^ = ^ g %3 Aidj 

In connection with the metric g, for vector fields u, v, 

*(u • v) = u b A (*v b ) 

so that consistently, 

(u,v) L 2 (Mj9) = (u b ,v b ) L 2 (Mi9) . (5.6) 

In a 2-dimensional Riemannin manifold, the divergence and curl of a vector field u € TM are 
then defined as scalar- valued functions^] on M, 

divergence, divu := — 5(u ) = *d * (u ) (5-7) 

curl, curlu := — * d(u t> ) (5-8) 

For a scalar field /, we define gradient and its tt/2 rotation as 

gradient, V/ := (dff (5.9) 

rotated gradient, V ± f := (£(*/))" = -{*dff (5.10) 

We also define the counterclockwise it/ 2 rotation operator - 1 acting on a vector field as 

u x := -(*u b )» (5.11) 

so that, consistantly, 

V ± / = (V/) ± , 

and 

div u = curl u -1 . 

It is then easy to use these definitions to verify the following properties, for scalar function / on 
a surface, 

curlV/ = divV^/ = (5.12) 

due to dd = and 55 = 0; and 

divV/ = curl V" 1 / = -5df (5.13) 



with — 5d being the classical surface Laplacian A (c.f. (5.4)). 



To this end, the vector-field version of Hodge decomposition (5.5) becomes 



u = V$ + V- L ^. (5.14) 



^Here and below, 0-forms are identified with scalar-valued functions. 
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We note that, by the virtue of (5.13), the decomposition satisfies 

divu = A<&, curlu = A*. 
We finally establish a version of the Green's identity and a version of the product rule on 



Riemannian manifolds. First, the duality relation (5.2) together with (5.4) and (5.6) implies 

(/, Ag) L2{M) = -(V/, Vg) L 2 {M>g y (5.15) 

Secondly, as a consequence of the product rule for differential d acting on wedge product, for scalar 
function z and vector field u, we have 



curl(zu ) = div (zu) = X7z ■ u + zdivu. 



(5.16) 



5.3. Local Expression in Terms of Spherical Coordinates for M = S 2 . Let (j) denote the 
logitude and 8 the colatitude of a point on a sphere. Let e^, eg denote the unit tangent vectors in 
the increasing directions of <j> and 6. Then, at point p that is away from the poles, 



do = e<j, 



namely, 



sin# 



: dj, and da 



form an orthonromal basis of TM„ 



(5.17) 



Therefore, the musical isomorphisms, in <f>, 6 coordinates, satisfy 



1 



sin 



"sin# 

In this context, the Hodge *-operator satisfies 



and (dg) = d6 form an orthonomal basis of T*M p 



for 1-forms, * (Aidtj) + A 2 d0) = -^—^dO - A 2 sin 

for 0-forms and 2-forms, * (Adcft A d6) 
The differential operators defined in (|5.7|) — (5.10) then become 



8111 ( 



A 



sm( 



*A = A sin Bdcj) A d6 



for vector field u =uie 9 j + u 2 eg 
1 



divu 



curlu 



sin I 
1 

sin( 



(dtpui + dg(u 2 sinO)) 
{dj>u 2 - dg(u\ sin 6)) 



and 



for scalar field / 
V/ = 



1 



8111 

V x f =dgfe, 
1 



+ 

1 

sin 6 



A/ 



sin 



(^/ + sin^(sin^/)) 
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6. Appendix B: Uniform Estimates Independent of e 



In this section, we use energy methods to prove local-in-time existence of classical solutions for 
the incompressible Euler equations independent of the Rossby number e. Rewrite the equation as 



in (2.6), 



1 



d t u + V^A" 1 curl (V u u) = ~£[u], 



where operator £, as in (2.5), is defined by 

£[u] := V ± A~ 1 curl (zu J 



(6.1) 



(6.2) 



The main challenge rises from the nontrivial geometry of S 2 : only a selective set of differential- 
integral operators on S 2 commute with each other. Although this is not a problem regarding 
well-posedness with fixed e, it causes difficulties in obtaining e-independent estimates. The fact 
that our £ has variable coefficients adds another layer of difficulties. In proving the following 
theorem, we will address these commutability issues specifically. 



Theorem 6.1. Consider the incompressible Euler equations (6.1), (6.2) on a rotating sphere S 2 
with div-free initial data uq. Given any integer k > 2, assume uq G H k {S 2 ). Then, there exists 
universal constants Co, To independent of e so that 

T 



||u(i, -)|| H fc < Colluolli^ for any t G 

Proof. For simplicity, we only prove the case when k is even. 
First, we show that for any u G L 2 (S 2 ), 



0. 



lUoHfffc 



/ 

Js 2 



u • £[u] = if divu = 0. 



(6.3) 



Indeed, by definition (6.2) and Hodge decompositon (2.1), (2.2), we have 

zu 1 - £[u] = (id - V^A" 1 curl )[zu ± ] = VA div [zu 1 ] 
which is curl-free and therefore L 2 -orthogonal to a div-free flow u, namely, 



Thus, 



/ u-(zu ± -C[u]) = Q. 
Js 2 

/ u-£[u] = / u- (zu 1 ) = 0. 
Js 2 Js 2 



Secondly, we show that A and £ commute. Indeed, for any incompressible flow u = V \P 
A£[u] = AV^A^ 1 ^ by (g§ 



V^AA"^^^ 
V ± A- 1 Ad 4 ,V 



by (5.4), (5.10) 
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The key step remaining is to show that A and 9^ commute. This can be done using the fact that 
spherical harmonics Y™ are eigenfunctions for both A and cL. More specifically, for any spherical 
harmonic 17™, 

= d^e^QT 1 (cos 6) = imYf 2 

and therefore 

Ad^™ = A(imY, m ) = -iml(l + 1)17" = d^AY™. 

Lastly, once the commutability of A and C are established as above, we easily obtain for even 
integer k > 0, 

f A k / 2 u • A fe / 2 £[u] = I A fc / 2 u • C[A k / 2 u] = 
Js 2 Js 2 

where the second equality is due to (|6.3[) and A fc ' 2 u also being incompressible (note: A and div 



commute). Now, take the I? inner product of A fc u with both sides of (6.1 ), knowing that the RHS 
should vanish, 

/ A fe u • d t u + A fc u • V ± A -1 curl (V u u) = 
Js 2 

and invoke the standard energy methods (e.g. [10J) to arrive at conclusion, which is clearly e- 
independent. □ 
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